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ABSTRACT 

1H 0707-495 is the most convincing example of a supermassive black hole with an X-ray 
spectrum being dominated by extremely smeared, relativistic reflection, with the additional 
requirement of strongly supersoler iron abundance. However, here we show that the iron fea¬ 
tures in its 2-10 keV spectrum are rather similar to the archetypal wind dominated source, 
PDS 456. We fit all the 2-10 keV spectra from 1H 0707-495 using the same wind model as 
used for PDS 456, but viewed at higher inclination so that the iron absorption line is broader 
but not so blueshifted. This gives a good overall fit to the data from 1H 0707-495, and an 
extrapolation of this model to higher energies also gives a good match to the NuSTAR data. 
Small remaining residuals indicate that the iron line emission is stronger than in PDS 456. 
This is consistent with the wider angle wind expected from a continuum driven wind from the 
super-Eddington mass accretion rate in 1H 0707-495, and/or the presence of residual reflec¬ 
tion from the underlying disk though the presence of the absorption line in the model removes 
the requirement for highly relativistic smearing, and highly supersoler iron abundance. We 
suggest that the spectrum of 1H 0707-495 is sculpted more by absorption in a wind than by 
extreme relativistic effects in strong gravity. 

Key words: black hole physics - radiative transfer - galaxies: active - galaxies: individual: 
1H0707-X-rays: galaxies. 


1 INTRODUCTION 

1H 0707-495 (hereafter 1H0707) is a Narrow Line Seyfert 1 
(NLS1) galaxy i.e. a low mass, high mass accretion rate (in terms of 
Eddington) AGN (Boroson 2002). It shows extreme dips in its X- 
ray lightcurve, during which the spectra have a steep drop around 
7 keV, associated with iron Ka (Boiler et al. 2002). These spec¬ 
tra are generally fit with ionized reflection, but the features are 
so strong and broad that this interpretation requires extreme con¬ 
ditions. The requirements are that the black hole spin is close to 
maximal, that the incident radiation is strongly focussed onto the 
inner edge of the disk whilst being suppressed in the direction of 
the observer, and that iron is overabundant by a factor of 7-20 
(Fabian et al. 2004, 2009, 2012; Zoghbi et al. 2010). The first two 
features can be explained together in a model where the dips are 
caused by an extremely compact X-ray source on the spin axis of 
the black hole approaching the event horizon (hereafter the lamp- 
post model). The resulting strong light bending focusses the in¬ 
trinsic continuum away from the observer (producing the drop in 


flux), whilst simultaneously strongly illuminating the very inner 
disk (Miniutti & Fabian 2004). 

However, the optical/UV continuum from 1H0707 im¬ 
plies that the black hole is accreting at a super-Eddington 
rate (Done & Jin 2016). Hence the inner disk is unlikely 
to be fiat, as assumed in the lamppost reflection models, 
and should launch a strong wind due to continuum radi¬ 
ation pressure (Ohsuga & Mineshige 2011; Jiang et al. 2014; 
Sjdowski & Narayan 2016; Hashizume et al. 2015). Thus this 
source, and other NLS Is with similarly very high mass accretion 
rates and similarly extreme X-ray spectra (e.g. IRAS 13224-3809: 
Ponti et al. 2010) may have spectral features which are affected by 
absorption and/or emission/scattering in a wind. Absorption has 
been persistently suggested as an alternative explanation for the 
features seen around iron in the NLS Is, though the spectra are com¬ 
plex and not easy to fit (Mkn 766: Miller et al. 2007; Turner et al. 
2007, MCG-6-30-15: Inoue & Matsumoto 2001; Gallo et al. 2004; 
Miller et al. 2008, 1H0707: Boiler et al. 2002; Miller et al. 2010; 
Mizumoto et al. 2014). 
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The recent discovery of blueshifted (v > 10000 km s _1 
i.e. 0.03c), narrow absorption lines from very highly ionized ma¬ 
terial (mostly He- and H-like Fe Ka) has focussed attention on 
winds from the inner disk. These Ultra-Fast Outflows (UFOs) are 
seen from a variety of nearby Seyfert galaxies (see e.g. the com¬ 
pilation by Tombesi et al. 2010). The fast velocities imply that 
the winds are launched from the inner disk since winds typically 
have terminal velocities similar to the escape velocity from their 
launch radius. However, the launch mechanism is not well un¬ 
derstood (e.g. Tombesi et al. 2012). Standard broad line Seyfert 
galaxies are well below their Eddington limit, so that they can¬ 
not power winds from continuum radiation pressure alone. More¬ 
over, they have strong X-ray emission which ionizes the disk ma¬ 
terial above where it has substantial UV opacity, suppressing a UV 
line driven disk wind (Proga et al. 2000; Proga & Kallman 2004; 
Higginbottom et al. 2014). Magnetic driving seems the only re¬ 
maining mechanism, but this depends strongly on the (unknown) 
magnetic field configuration so is not yet predictive (Proga 2003). 


Despite this general lack of understanding of the origin of 
the winds, the fastest, v > 0.2c, and most powerful UFOs are 
seen in luminous quasars such as PDS 456 (Reeves et al. 2009; 
Nardini et al. 2015) and APM 08279+5255, a high redshift quasar 
which is gravitationally lensed (Chartas et al. 2002). These AGN 
are not standard broad line Seyfert galaxies. Instead, since both 
quasars are around the Eddington limit, they could power a con¬ 
tinuum driven wind. Also, black holes in both quasars are very 
high mass (~ 10 9 — 10 10 Mq) so that their accretion disk spec¬ 
tra should peak in the UV. PDS 456 is also clearly intrinsically 
X-ray weak (Hagino et al. 2015), and APM 08279+5255 may also 
be similar (Hagino et al, in preparation). A more favorable set of 
circumstances for UV line driving is hard to imagine. Nonetheless, 
the extremely high ionization of the UFO in PDS 456 means that 
the observed wind material has no UV opacity, so UV line driving 
must take place out of the line of sight if this mechanism is impor¬ 
tant (Hagino et al. 2015). 


Whatever the launch mechanism, PDS 456 clearly has features 
at Fe Ka which are dominated by a wind from the inner disk rather 
than extreme reflection from the inner disk (Reeves et al. 2003, 
2009; Nardini et al. 2015). The best wind models so far use Monte- 
Carlo techniques to track the complex, geometry and velocity de¬ 
pendent, processes in the wind including absorption, emission, con¬ 
tinuum scattering and resonance line scattering (Sim et al. 2008, 
2010; Hagino et al. 2015). While these models have been used to 
fit Mkn 766 (Simetal. 2008), PG 1211 + 143 (Simetal. 2010), 
PDS 456 (Hagino et al. 2015), and six other ‘bare’ Seyfert 1 nu¬ 
clei (Tatum et al. 2012), they have not yet been applied to 1H0707, 
the object with the strongest and broadest iron features, and the one 
where the reflection models require the most extreme conditions 
(Fabian et al. 2009, 2012). Previous work showed that simple wind 
models, where the iron features were described using a single P 
Cygni profile, can fit the deep drop at ~ 7 keV seen in one ob¬ 
servation of this object (Done et al. 2007), but here we use the full 
Monte-Carlo wind code, MONACO, of Hagino et al. (2015) to see if 
this can adequately fit all the multiple observations of 1H0707. We 
show that inner disk wind models can indeed give a good overall 
description of the 2-10 keV spectra, and that extrapolating these to 
higher energies can also match the NuSTAR data from this source. 


Table 1. XMM-Newton and Suzaku observations of 1H0707 


Name 

Obs ID 

Start Date 

Net exposure (ks) 1 

XMM-Newton 

Obsl 

0110890201 

2000-10-21 

37.8 

Obs2 

0148010301 

2002-10-13 

68.1 

Obs3 

0506200301 

2007-05-14 

35.8 

Obs4 

0506200201 

2007-05-16 

26.9 

Obs5 

0506200501 

2007-06-20 

32.6 

Obs6 

0506200401 

2007-07-06 

14.7 

Obs7 

0511580101 

2008-01-29 

99.6 

Obs8 

0511580201 

2008-01-31 

66.4 

Obs9 

0511580301 

2008-02-02 

59.8 

ObslO 

0511580401 

2008-02-04 

66.6 

Obsl 1 

0653510301 

2010-09-13 

103.7 

Obsl 2 

0653510401 

2010-09-15 

102.1 

Obsl 3 

0653510501 

2010-09-17 

95.8 

Obsl 4 

0653510601 

2010-09-19 

97.7 

Obsl 5 

0554710801 

2011-01-12 

64.5 

Suzaku 

SuzakuObs 

700008010 

2005-12-03 

97.9/100.4/97.2/97.8 


1 Net exposure time of pn for XMM-Newton and XIS0/XIS1/XIS2/XIS3 
for Suzaku , respectively. 



Figure 1. All 2—10 keV spectra observed by XMM-Newton /pn and 
Suzaku! FI detectors, unfolded against a F = 2 power law. 


2 OBSERVATIONS OF 1H0707 AND COMPARISON TO 
THE ARCHETYPICAL WIND SOURCE PDS 456 

1H0707 was observed by Suzaku and XMM-Newton for many times 
as listed in Tab. 1 . We reduce both XMM-Newton pn and Suzaku 
XIS data with standard screening conditions: PATTERN < 4 
events for XMM-Newton/ pn data and grade 0, 2, 3, 4 and 6 events 
for Suzaku/XlS data were used. Bad time intervals were also ex¬ 
cluded. XMM-Newton data in time intervals when background 
rates of PATTERN = 0 events at energy > 10 keV exceeded 
0.4 counts s- 1 were removed. Suzaku data within 436 s of pas¬ 
sage through the South Atlantic Anomaly (SAA), and within an 
Earth elevation angle (ELV) < 5° and Earth day-time elevation 
angles (DYE_ELV) < 20° were excluded. Spectra were extracted 
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Figure 2. Comparison between the brightest/faintest spectra of PDS 456 
(black/grey: ObsID 707035030/701056010) and the brightest/faintest spec¬ 
tra of 1H0707 (cyan/blue: Obsl2/Obsl5). The 2-10 keV spectra show a 
similar range in 2-5 keV spectral slope, and similarly deep features at Fe 
K a, though these start at a lower energy in 1H0707, and are broader than in 
PDS 456. 


from circular regions of 64” and 2/9 diameter for XMM-Newton 
and Suzaku, respectively. Background spectra of XMM-Newton 
data were extracted from circular regions of the same diameter in 
the same chip as the source regions, while background spectra of 
Suzaku data were extracted from annular region from 7/0 to 14/0 
diameter. The observed data show a large variability in continuum 
spectra as seen in Fig. 1. There is an obvious variability in the 2- 
6 keV continuum shape and in the strength of the iron K features 
around 7 keV. 

Figure 2 shows brightest and faintest 2-12 keV (rest frame) 
spectra seen in PDS 456, the archetypal inner disk wind source, 
(black/grey: ObsID 707035030/701056010), together with the cor¬ 
responding brightest/faintest observations of 1H0707 (cyan/blue: 
Obsl2/Obsl5). The spectra are similarly harder when fainter, and 
show a similarly deep drop at iron K. In PDS 456, these features 
are all generally associated with absorption, with the deep drop 
at iron K interpreted as Fe Ka absorption lines from highly ion¬ 
ized material in the wind, while the variability at lower energies 
can be described by complex absorption from lower ionization ma¬ 
terial (Reeves et al. 2003, 2009; Nardini et al. 2015; Hagino et al. 
2015). Conversely, the spectral features in 1H0707 are generally as¬ 
sociated with extreme relativistic reflection, especially in the hard¬ 
est/dimmest spectra (Fabian et al. 2009, 2012). Nonetheless, the 2- 
10 keV spectra of 1H0707 and PDS 456 are quite similar, both in 
range of continuum shapes and in the shape of the drop around Fe 
Ka. 

However, there are also some more subtle differences. The 
sharp drop in 1H0707 is at ~ 7.0-7.5 keV, whereas it is at 
8.5-9.0 keV in PDS 456, and the recovery of the continuum after 
the absorption lines at higher energies is more evident in PDS 456 
than in 1H0707. Thus if these features are from a wind in 1H0707, 
it shares many similarities with the wind in PDS 456 but cannot 
have exactly the same parameters. 


3 THE MONACO SIMULATIONS FOR PDS 456 

The wind in PDS 456 was modeled by using the MONACO Monte- 
Carlo code. This wind model and the calculation scheme are fully 
described in Hagino et al. (2015), and the general framework de¬ 
sign of the code is described in Odaka et al. (2011). To summarize 
its main features, it follows Sim et al. (2010) and assumes a bicon- 
ical wind geometry, where the wind is launched from radii R m i n to 
Rm ax along streamlines which diverge from a focal point at a dis¬ 
tance d below the black hole. The radial velocity of material along 
each streamline of length l from its launch point on the disk is 

Vr{l) = Vo + (Uoo - Vo) ( 1 - Rmn " \ , (1) 

\ 4lmin it/ 

while the azimuthal velocity is assumed to be Keplarian at the 
launch point, and then conserves angular momentum as the wind 
expands. We assume vo is negligible. There can also be a turbulent 
velocity, vt. 

The total wind mass loss rate M w i n d is given as 

ATwind = 1.23m p n(r)v r (r)4nD 2 -f, (2) 

47T 

which sets the density n(r ) as a function of radius in the wind. 
Here, 1.23m p , D and fi are an ion mass, the distance from the fo¬ 
cal point and a solid angle of the wind, respectively. The bicone 
is split into 100 radial shells, and the ion populations are calcu¬ 
lated in each shell by using the XSTAR photoionization code on the 
assumption that the central source is a power law spectrum with 
photon index T and ionizing luminosity L x (adjusted for special 
relativistic dimming from the radial outflow velocity). The result¬ 
ing H- and He-like ion densities of each element are put in the 
Monte-Carlo radiation transfer simulator MONACO, which uses the 
GEANT4 toolkit library (Allison et al. 2006) for photon tracking, 
but extended to include a full treatment of photon processes related 
to H- and He-like ions (Watanabe et al. 2006). 

Thus there are 9 major free parameters. R, n , n , -Rmax and d de¬ 
termine the geometry. Vao, Vt and /3 define the radial velocity struc¬ 
ture. These geometry and velocity structure together with M„ i n d 
determine the density structure, which then sets the ionization state 
given L x and T. 

This wind model is self-similar in ionization structure and col¬ 
umn density for systems at different mass but the same Eddington 
ratio m = Tboi/T/Edd he. the same geometry, velocity structure 
and spectral index. The ionization state £ = L x /nD 2 can be writ¬ 
ten as 

£ = 1.23m p u r DL x /M W md, (3) 

where equation 2 is used. Since the ionizing luminosity can be writ¬ 
ten as L x = f^thr/MEddC 2 where / x = L x /Lboi, then the ioniza¬ 
tion parameter reads 

£ = 1.23m p v r D/ x mr7c 2 /(M win d/M E dd). (4) 

Thus the ionization state of the wind is determined by 
ATwind/ATEdd, and not by black hole mass directly. 

The observed spectral properties strongly depend on the in¬ 
clination angle. The energy at which the line absorption is seen 
depends on both the terminal wind speed and the inclination an¬ 
gle at which we see the wind, whereas the width of the absorption 
line depends on the spread of velocities along the line of sight. This 
spread in turn depends on the terminal velocity of the wind, how 
much of the acceleration region is along the line of sight, and the 
projected angle between the wind streamlines and the line of sight. 
Along the bicone, a fast wind acceleration law means that most of 
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the wind is at its terminal velocity. Thus the line width is fairly 
small and the blueshift indicates the true wind velocity. At higher 
inclination angles, the line of sight cuts across the acceleration re¬ 
gion so the absorption line is wider and the total blueshift is not so 
large since the line of sight includes much lower velocity material. 
Thus the same wind seen at different line of sight can have very 
different properties in terms of the measured width and blueshift of 
the absorption line. 

We illustrate this inclination dependence in Fig. 3 by showing 
the MONACO wind model for PDS 456 at a larger range of incli¬ 
nation angles than in Hagino et al. (2015). This wind model has 
_R min = d = 20-Rg and R m ax = 1.5f? m in so that the wind fills 
a bicone between 0 m i n = 45° to # max = 56°.3, i.e. the wind fills 
a solid angle fl/47r = 0.15. We fix the velocity law at /? = 1 
and assume vq = Vt — 10, 000 km/s. The terminal velocity is 
set at Voo = —0.3c, which is consistent with the escape velocity 
from -Rmin- The mass outflow rate is M w ind/A^Edd = 0.3, and the 
ionizing spectrum with T x /I/Edd = 10 44 /10 4 ' and T = 2.5 is 
assumed. It is clear that inclination angles larger than 0 max gives a 
line which is less blueshifted, but broader, as required to fit 1H0707. 
The width is accentuated because the larger column density of the 
wind means that the ionized iron edge is also important in absorb¬ 
ing the spectrum blueward of the absorption line. These high incli¬ 
nation MONACO model spectra for PDS 456 give iron K features 
which are quite similar to those required for 1H0707. 

4 FITTING THE RANGE OF 2-10 KEV DATA FROM 
1H0707 

4.1 Tailoring the wind model to 1H0707 

We first see whether the MONACO model calculations for the 
archetypal wind source PDS 456 can fit the features seen in 1H0707 
by simply viewing the wind at a larger inclination angle. We 
fit to Obsl2, which is the steepest spectrum seen from 1H0707. 
We assume that this steepest spectrum has negligible absorption 
from lower ionization species. We use the MONACO results from 
PDS 456, as coded into a multiplicative model by Hagino et al. 
(2015), on a power law continuum. Free parameters are the incli¬ 
nation and the redshift. By allowing the redshift to be free, we are 
able to fit for a slightly different wind velocity. This fitting gives 
Xv 2 = 133/87 for a power law index of 2.65 ± 0.05 (see Fig. 4a 
and Table 2). 

We compare this model with the standard extreme reflection 
interpretation. We use the ionized reflection ATABLE models RE- 
FLIONX (Ross & Fabian 2005), and convolve these with KDBLUR, 
allowing the emissivity index, inner disk radius, inclination angle, 
iron abundance and ionization parameter to be free parameters. 
This gives \„ 2 = 115/83 for T = 3.07 ± 0.15, a somewhat bet¬ 
ter fit, but with more free parameters (Fig. 4b). Inspection of the 
residuals shows that the wind model underpredicts the emission 
line flux, whereas the reflection model underpredicts the extent of 
the drop at ~ 7.5 keV. This suggests that the drop is better matched 
by the wind, but there is more line emission than predicted in the 
PDS 456 models. 

This discrepancy could be produced by the wind itself. A 
wider angle wind will intercept more of the source luminosity, and 
have stronger emission/reflection/scattering features at iron K for a 
given absorption line strength. We will investigate such wide angle 
winds in a subsequent paper. Alternatively (or additionally), there 
can be lower ionization species in the wind which also add to the 


line emission. We approximate both of these physical pictures by 
adding a phenomenological broad gaussian line to the wind model, 
and find a better fit than either reflection or the PDS 456 wind alone 
(X 2 = 107/84: Fig. 4c). 

Another possible origin for the excess iron line emission is 
that there is residual reflection from the underlying disk. We add a 
phenomenological blurred reflection component as above, allowing 
the amount of disk reflection, emissivity index, inner disk radius, 
inclination angle, iron abundance and ionization parameter to be 
free. This is a better fit than the gaussian line (x 2 = 92.0/82), 
mostly because the reflection model also includes hydrogen-like 
sulphur line emission at 2.8 keV which is clearly present in the data. 
However, the best fit reflection parameters are now more extreme 
than those for the fit without a wind, with strongly centrally peaked 
emissivity and fairly small inner radius (Table 2). 

Such strongly centrally peaked emissivity is characteristic of 
the lamppost model, but this produces an illumination emissivity 
profile which is more complex than a single power law. Hence 
we replace KDBLUR with the fully relativistic lamppost model for 
the emissivity (RELCONV_LP; Dauser et al. 2013). Fixing the black 
hole spin to maximal and allowing the inner radius to be free (as 
in KDBLUR) gives a slightly worse fit, with xt — 94.4/82 for 
similarly extreme parameters (inner radius of ~ 3 r g and height 
of < 2 r g : Fig. 5a). However, this fit is very poorly constrained 
as shown in Fig. 5b. Only the most extreme solution at small Ri n 
and small h can be excluded with more than 99% confidence level, 
while the best fit has very similar fit statistic for much less extreme 
fits (Ri n ~ 30 r g for height which is unconstrained). We illustrate a 
statistically equivalent less extreme fit in Fig. 4d by using KDBLUR 
with fixed emissivity of 3 and find a much larger inner radius of 
~ 36 r g (\t = 95.1/83: Table 2). This result shows directly that 
with wind absorption, the lamppost reflection is not required to be 
extreme in either relativistic smearing nor in iron abundance. 

Regardless of how the additional emission component is mod¬ 
eled, our wind model developed for the highly blueshifted absorp¬ 
tion lines in PDS 456 is better able to produce the observed sharp 
drop at ~ 7 keV in 1H0707 than relativistic reflection alone. 

4.2 Application to XMM-Newton and Suzaku data 

We re-run MONACO with the parameters around the values obtained 
in the previous section but now put in explicitly the mass of 5 x 
10 6 Mq, an observed X-ray spectrum with T = 2.6 and L x = 
2 x 10 42 erg s^ 1 , and reduce the wind velocity to 0.2c but keep the 
same velocity law. 

The result from the previous section suggests that there is 
rather more line emission compared to absorption than is pre¬ 
dicted for the PDS 456 wind geometry. However, a larger solid 
angle of the wind would necessitate more resolution in 9 and a 
change in the way of calculating the ionization state. Hence we 
keep the same wind geometry for the present, and run simulations 
for M w i n d/MEdd = 0.13, 0.20, 0.26. Refitting these to Obsl2 
gives x 2 = 139.7,132 and 133 for 89 degrees of freedom, re¬ 
spectively. Hence we pick M w i n d/MEdd = 0.2, and apply this 
model to all the observations of 1H0707. Similarly to PDS 456, we 
assume that the intrinsic power law stayed constant in spectral in¬ 
dex, and that the major change in spectral hardness is produced by 
a changing absorption from lower ionization species. To describe 
this absorption, we used a very simple model of partially ionized 
absorption which can partially cover the source (ZXIPCF in XSPEC). 
This gives an adequate fit to all the PDS 456 spectra, including the 
hardest (black line in Fig. 2) which bears some similarities to the 
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45-47 deg 51-53 deg 57-59 deg 








Figure 3. MONACO wind model for PDS 456 seen at different inclination angles. The direct and reprocessed components are plotted in red and blue, respec¬ 
tively. The sums of these components are plotted in black. Angles from 45 — 47° intercept only the top of the bicone. Almost all the wind is at the terminal 
velocity of 0.3c so the lines are strongly blueshifted, but fairly narrow. Larger angles cut through more of the wind, so the lines are deeper, and intersect 
progressively more of the acceleration zone, so the lines are broader and blend into each other. Angles larger than 56.6° now cut across the bicone, so the lines 
are less blueshifted, but include even more of the acceleration zone, so the lines are broader, and now the wind is close to optically thick so the edge is also 
important. Thus the line absorption shifts from two obviously narrow lines which are strongly blueshifted, to a single, very broad absorption line which is not 
strongly blueshifted. 


Table 2. Fitting parameters for Obsl2 with MONACO wind models and the extremely blurred reflection model. 





(b) KDBLUR 

(c) MONACO wind 

(d) MONACO wind 



(a) MONACO wind 

*REFLIONX 

+ ZGAUS 

+ KDBLUR* 

REFLIONX 






Emissivity:free 

Emissivity:3 

MONACO wind 

Velocity (c) 

0-19±o!ol 

— 

> 0.17 

p 

00 
l + 
0 9 
b b 

to to 

p 

00 

1 + 
P P 
b b 

to 


dincl (°) 

52.6± 2 ; 7 

— 

62 4”*” 2,0 

10 9 

51.4t 2 ;? 

52.4^3 ® 

Powerlaw 

r 

z.uO—004 

W 

b 

-4 

1 + 
o o 

to 00 

9 oq+0.09 
z.oo_oio 

2.89±g;ig 

9 +0+9 

KDBLUR 

Index 

_ 

Q C 7 +O .66 
0.01 _0.35 

_ 

> 4.39 

3.0 1 


-^in g) 

— 

< 10.85 

— 

q q-|-0.6 

o.O —06 

36i 7 r 


dincl (°) 

— 

48.901 3 ; 3 ® 

— 

tied to wind 

tied to wind 

REFLIONX 

Fe abundance 

— 

> 6.75 

— 

q ni +4.79 
o.Ul-i es 

1 cifi+5-72 
-L-OO-i 36 


? 

— 


— 

list *® 7 

6 QQ + 534 

oyy -650 

Gaussian 

Line energy 

— 

— 

6.68l“:® 3 

— 

— 


Sigma 

— 

— 

> 1.18 

— 

— 

Fit statistics 

X 2 /dof 

133.4/87 

115.5/83 

107.4/84 

92.0/82 

95.1/83 


Null prob. 

1.0 x 10“ 3 

1.1 x 10~ 2 

4.4 x 10~ 2 

0.21 

0.17 


1 Parameters are fixed. 


hardest spectra in 1H0707, which are generally interpreted as re¬ 
flection dominated (Fabian et al. 2012). 

All the spectra (black) and models (red) are plotted in Fig. 6. 
We also show the background (blue) for each spectrum, so that the 
signal-to-noise at high energies can be directly assessed. The pa¬ 
rameters for each fit are listed in Tab. 3. For Obs8 and Obsl2, the 
absorption component is not included because the significance of 


adding ZXIPCF does not exceed the 99% confidence level of an F- 
test. 

As shown in Fig. 6, our MONACO models roughly reproduce 
the sharp drop around 7.1-7.5 keV in all the observed spectra. The 
best fit values of inclination angle ranges from ~ 58° to 75°. These 
are all higher than the inclinations inferred for PDS 456, and there 
is a broader spread in derived inclination angle. 
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Figure 4. XMM-Newton /pn spectra of Obsl2 fitted by (a) WlND*POWERLAW (top-left), (b) POWERLAW+KDBLUR*REFLIONX (top-right), (c) 
wind(powerlaw+ZGAUS) (bottom-left) and (d) wind(powerlaw+kdblur*reflionx) with emissivity fixed at 3 (bottom-right). The spectra and models 
are plotted in black and red, respectively. The lower panels show the residuals in units of x- 



Figure 5. XMM-Newton /pn spectra of Obsl2 fitted by the wind model with lamppost disk reflection (WIND(P0WERLAW+RELC0NV_LP*REFLI0NX), left 
panel) and contour plot of the inner disk radius and the height of the X-ray compact source in the lamppost geometry (right panel). The red, green and blue 
lines in the contour plot correspond to confidence levels of 68, 90 and 99 percent, respectively. The color map shows the value of % 2 for each parameter grid. 
The cyan cross is plotted at the best fit, where h ~ 1 r g and Ri n ~ 3r g . The model spectrum shown in the left panel corresponds to this extreme solution. 
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Figure 6. Observed spectra and MONACO models with v = 0.2c and M w i n d/MEdd = 0.2 for all the observations. The spectra and models are plotted in 
black and red, respectively. The background spectra are also overplotted in blue. The lower panels show the residuals in units of x- 
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Table 3. MONACO fit with v = 0.2c and M w ind/-^Edd = 0.2 to all the observations of 1H0707. 


Name 

Continuum absorption 


MONACO wind 

Fit statistics 


A h (10 22 cm -2 ) 

log? 

/cov 

^incl 

X 3 /dof 

Null probability 

Obsl 

28l® 3 

1 90"*”°' 91 

-L.tnj—o 63 

0 «5+ 0,12 
U.OO_ 0 03 

74.61 3 '® 

34.5/40 

0.73 

Obs2 

66 1« 

1 90”*" 0 ' 28 

0-45lg;« 

64.ll 2 ; 3 

85.2/79 

0.30 

Obs3 

isir 

< 1.63 

0-681^5 

64.il 3 ;® 

52.7/55 

0.56 

Obs4 

o 9+18.2 

° -6.3 

< 2.28 

0-80l 3 ;Jg 

57.9 1 

15.1/21 

0.82 

Obs5 

9q+72 
zo —16 

< 2.14 

O oq+ 0-44 
U-^o_o 09 

68.0+g 

102.6/69 

0.005 

Obs6 

o 4 +90.6 
8 - 4 -4.5 

< 2.02 

0 42+ 0 - 29 
u *^ z —0.14 

63.4lf 2 6 7 

28.2/34 

0.75 

Obs7 

iSltw 

9 71 

* —0.56 

0 71+0.06 
u>1 ± -0.25 

68.3l^ 

115.1/87 

0.02 

Obs8 

— 

— 

— 

64.0+ 3 ® 

129.12/84 

0.001 

Obs9 

65t« 

< 2.26 

0.42l°;« 

58.3l 3 g 

88.6/79 

0.22 

ObslO 

124 ±H 9 

9 71 “bO- J-9 
‘■ L —1.90 

0-58l°;l| 

68 .21 22 

87.4/72 

0.10 

Obsl 1 

121+32 

1 Z1_ 73 

2 76 +0 - 21 

' D —0.46 

0 60”*”°' 20 
U.OU-o 24 

63 1~*~ 3,3 
90.1—2 e 

100.3/81 

0.07 

Obsl 2 

— 

— 

— 

60.6+g 

132.0/89 

0.002 

Obsl 3 

H7+ 28 

11 '-63 

O 7 o + 0.15 

1 0.31 

0-63l^ 3 

63.21 2 ! 

154.9/78 

5 x 10 -7 

Obsl4 

181±g 03 

o OQ + 0.10 

z,.°^_o 09 

0-82lg;°s 

63.ll 33 

75.3/77 

0.53 

Obsl 5 

23 If 

1 61”*” 1 ' 02 
±0± -0.78 

o.93lg;g* 

70.91 ll 

43.1/51 

0.78 

SuzakuObs 

109125 

o 7 9+0.07 
z -‘ z -1.47 

o.89l“;“ 

67.o!g4 

97.9/86 

0.18 


1 No constraints on this parameter were obtained. 
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Figure 7. NuSTAR data compared with the best-fit wind model of Obsl5. 
The spectra of XMM-Newton, NuSTAR are plotted in grey and black, while 
the model lines for XMM-Newton and NuSTAR are plotted in magenta and 
red, respectively. The lower panels show the residuals in units of 


4.3 Application to NuSTAR data 

A key breakthrough in AGN spectral studies has come from NuS¬ 
TAR data which extend the energy range of the observed spectra 
beyond 10 keV and have much better signal-to-noise above 7 keV 
than the Suzaku or XMM-Newton data (see Fig. 6). 1H0707 has also 
been observed by NuSTAR, and the resulting spectra can be fit by 
the extreme relativistic reflection models (Kara et al. 2015). 

We extract the NuSTAR data, and follow Kara et al. (2015) in 
selecting the second NuSTAR dataset, which is a good match to 
XMM-Newton Obsl5. Fig. 7 shows the model fit to Obsl5 extrap¬ 
olated to 30 keV. We note that around half the drop at ~ 7 keV is 
from the wind model alone, while the other half is from the com¬ 
plex lower ionization absorption. There are no additional free pa¬ 
rameters, but the model gives a good fit to the higher energy data. 


5 DISCUSSION 

5.1 Effects of cool clumps 

We model the wind using the MONACO Monte-Carlo code, which 
tracks emission, absorption and scattering in a continuous wind ge¬ 
ometry (Hagino et al. 2015). In this continuous wind model, all the 
atoms lighter than iron are almost fully ionized, so that it cannot 
produce the strong spectral variability seen below the iron line re¬ 
gion. Instead, in the archetypal wind source PDS 456, such variabil¬ 
ity is assumed to be from lower ionization species. In PDS 456 this 
additional spectral variability can be approximately modeled us¬ 
ing partially ionized material which can partially cover the source. 
We show that this same combination of partial covering by lower 
ionization material together with the highly ionized wind gives an 
acceptable fit to the 1H0707 spectra, including the higher energy 
data from NuSTAR. 

The lower ionization material in the wind is likely to be 
clumped due to the ionization instability for X-ray illuminated ma¬ 
terial in pressure balance (Kroliketal. 1981). These clumps are 
cooler and less ionized than the hot phase of the wind, but with 
lower filling factor. The dramatic dips which are characteristic of 
the light curves of complex NLS1 like 1H0707 can then be in¬ 
terpreted as occultations by these cool clumps. We note that time 
dependent, clumpy absorption is typical of both UV line driven 
disk winds (Proga & Kallman 2004), and continuum driven winds 
(Takeuchi et al. 2014). 

The partially ionized clumps which partially cover the source 
imprint some spectral features as well as curving the continuum. 
Fig. 8 shows the data and model in the entire energy band for the 
highest (top panel) and lowest (bottom panel) flux states. The data 
is plotted in black and the total model spectrum is shown by the red 
line. The model spectrum is separated into absorbed (blue), unab¬ 
sorbed (green) and a soft X-ray excess component (magenta, see 
below). The low ionization absorption imprints strong atomic fea¬ 
tures from oxygen and iron-L below 2 keV in the lowest intensity 
state, but these are diluted by the X-rays which are not absorbed, 
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so that the total spectrum is almost featureless in the 0.5-2 keV 
bandpass. Hence it predicts no observable features in the RGS. 

The unabsorbed component does not dilute the iron-K features 
because the absorbed component is dominant at higher energies. 
The effect of absorption at iron K from the cool clumps alone is 
shown by the pink line in the inset window of Fig. 8, while the 
red line in the inset shows the total (cool plus hot phase) absorp¬ 
tion. The pink line shows that the cool clumps do imprint iron Ka 
absorption lines around 6.45 keV rest frame from Fe XVIII-XX as 
well as the much stronger K/3 lines around 7 keV rest frame. The 
K/3 lines blend with the hot phase absorption, while the Ka lines 
have an equivalent width of only a few eV, which is not detectable 
even with a 300 ks observation with Hitomi for a source as weak as 
1H0707. Hence the model is consistent with current and near future 
limits on spectral features in the data. 

A soft excess component (magenta) is added in the model 
spectrum in Fig. 8 in order to roughly describe the spectrum at 
low energies. We use a thermal Comptonization model COMPTT, 
allowing only the normalization to be free while fixing the shape 
to parameters typically seen in other AGN, namely a seed photon 
temperature of 0.05 keV, coronal temperature of kT = 0.2 keV and 
optical depth at r = 15 (Done et al. 2012). The plots show that the 
observed spectrum is roughly reproduced by this model, but there is 
an excess around 0.8 keV and a deficit between ~ 1-2 keV. These 
differences can be interpreted as broad emission and absorption (P 
Cygni profile) from oxygen K and iron L-shell transitions, which 
are naturally expected in the ionized fast wind, but are not included 
in our current Monte-Carlo code. Hence these features in the data 
are not well described by the model used here, which is tailored 
instead to iron K. Nonetheless, we do expect that similar features 
would be produced by including emission and absorption from cool 
clumps in our wind model. 

Our model is currently incomplete as these clumps should 
also reflect/emit as well as absorb, adding to the reflection/emission 
from the highly ionized phase of the wind which is included in our 
MONACO simulation. This reflection/emission from the composite 
wind could dominate during the dips, where the cool clumps prob¬ 
ably completely cover the source itself. We also envisage the cool 
clumps to be entrained in the wind, so the atomic features should 
also be blueshifted/broadened by the wind velocity structure. This 
means that the cool phase of the wind also contribute to the total 
wind kinetic luminosity which has so far been neglected. 

The clumpy structure can also play an important role for 
a characteristic time lags observed in 1H0707. The observations 
show a complex pattern of lags between soft and harder energies, 
with the soft band leading for slow variability, and lagging for 
fast variability (Fabian et al. 2009; Zoghbi et al. 2010). In the rela¬ 
tivistic reflection picture, this can be explained by a partially ion¬ 
ized disk, where reflection is weaker at 2-4 keV than at lower or 
higher energies. Thus, reflection dominates in the soft band, while 
the intrinsic continuum dominates in the intermediate band at a 
few keV. Reflection follows the continuum, but with a lag from 
the light travel time from the source to the reflector, giving the 
soft (conventionally referred to as a ‘negative’) lag at high fre¬ 
quencies. The positive lag at lower frequencies is probably from 
propagation of fluctuations through the accretion flow. On the other 
hand, a full spectral-timing model including both the positive and 
negative lags for PG 1244+026, another AGN with very similar 
mass and mass accretion rate to 1H0707, gives the observed lag of 
200 s seen in PG 1244+026 from reprocessing by material extend¬ 
ing from 6-12i? g (assuming a mass of 10'Mg: Gardner & Done 
2014). In this interpretation, this lag is due to the illuminating 



0.5 1 2 5 10 20 30 

Energy (keV) 


Figure 8. The observed data and model of the highest flux state (Top: 
Obsl2) and the lowest flux state (Bottom: Obsl5 and the NuSTAR data) 
plotted over a wide energy range. The absorbed spectrum is shown as the 
blue line, the unabsorbed in green, the total in red, and the observed data in 
black. A soft X-ray excess component with typical parameters is also plot¬ 
ted in magenta to show the spectrum at lower energies below 2 keV. There 
are no spectral features predicted below 2 keV due to the dilution by the 
unabsorbed (or reflected/scattered) flux (green). There are absorption lines 
from Fe XVIII—XX produced by cool clumps (pink line in a small window 
inside bottom panel) but the Ka lines at ~ 6.45 keV are not detectable even 
with 250 ks of Hitomi, and the higher energy lines merge with the hot phase 
of the wind (red line in a small window). 


hard X-ray flux which is not reflected but thermalized and re¬ 
emitted as quasi-blackbody radiation. Adding multiple occultations 
by clumps onto the same full spectral-timing model reduces the 
predicted lags to 50 s, similar to the 30 s observed in 1H0707 
but without requiring a small reflector distance (Gardner & Done 
2015). Thus cool clumps may be able to explain both the time aver¬ 
aged spectrum and the extremely short lags without requiring light¬ 
bending from a small corona close to the event horizon of an ex¬ 
treme spin black hole, though it remains to be seen whether they 
can fit the lag-energy spectra as well as the lag frequency spec¬ 
tra (Gardner & Done 2015). However, we note that the lamppost 
model itself also has difficulty fitting the details of the lag-energy 
spectra (Wilkins et al. 2016). Their alternative geometry of an ex¬ 
tended corona may give a better fit to the lag-energy spectra, but it 
seems unlikely to be able to simultaneously explain the deep dips 
in the light curves which can be explained in the lamppost model 
by the very small source approaching the horizon. 

5.2 Wind geometry 

Our specific wind model is probably not unique in terms of wind 
geometry and velocity structure. In particular, we used the same 
wind solid angle of fl/47r = 0.15 in 1H0707 as for PDS 456. 
This is likely to be appropriate for PDS 456 if there is indeed some 
component of UV line driving to the wind. It is because numerical 
studies show that this driving mechanism results in time variable 
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but fairly narrow wind streamlines. However, continuum radiation 
driving is much more likely in 1H0707, in which case the numer¬ 
ical simulations predict a wider opening angle wind. Future work 
should simulate winds with a larger opening angle, which would 
require correspondingly larger mass loss rates in order to keep the 
same absorption column density. 

The wider opening angle wind should also result in stronger 
line emission, so that this may give a better fit to some of the re¬ 
maining residuals. However, we caution that all these spectra are 
co-added, integrating over dramatic short timescale variability so 
steady state models may not be appropriate. 


6 CONCLUSIONS 

We can successfully reproduce the range of spectra seen from 
1H0707 with absorption in a wind from the inner disk. The strong 
spectral drop is produced by our line of sight cutting across the 
acceleration region where the wind is launched. In this region, we 
see the absorption line over a wide range of velocities, making a 
very broad feature. Since the column density in the wind is large, 
the broad absorption line merges into the absorption edge, which 
further depresses the spectrum at higher energies. 

Our wind model probably inappropriate to this source as dis¬ 
cussed in section 5. The time variability and the low ionization ab¬ 
sorption in the observed spectra seems to be consistent with an ex¬ 
istence of cool clumps in the wind. This clumps are naturally ex¬ 
pected due to the instability of the wind. Also, a larger opening an¬ 
gle of the wind would probably give better fits to the observed data. 
This implies the continuum radiation driving mechanism, which is 
expected in the super Eddington accretion. 

Nonetheless, the wind model presented here can fit the over¬ 
all 2-30 keV spectra from 1H0707. Unlike the lamppost models, it 
does not require that the black hole has extreme spin, nor does it re¬ 
quire that we have a clean view of a flat disk, nor does it require that 
iron is 7-20 times overabundant. We suggest that the extreme iron 
features in 1H0707 arise more from absorption/scattering/emission 
from an inner disk wind than from extreme gravity alone. 
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